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Two-Body Dirac equations of constraint dynamics provide a covariant framework to investigate 
the problem of highly relativistic quarks in meson bound states. This formalism eliminates automat- 
ically the problems of relative time and energy, leading to a covariant three dimensional formalism 
with the same number of degrees of freedom as appears in the corresponding nonrelativistic prob- 
lem. It provides bound state wave equations with the simplicity of the nonrelativistic Schrodinger 
equation. Unlike other three-dimensional truncations of the Bethe-Salpeter equation, this covariant 
formalism has been thoroughly tested in nonperturbatives contexts in QED, QCD, and nucleon- 
nucleon scattering. Here we continue the important studies of this formalism by extending a method 
developed earlier for positronium decay into two photons to tests on the sixteen component quarko- 
nium wave function solutions obtained in meson spectroscopy. We examine positronium decay and 
then the two-gamma quarkonium decays of rjc,r/' c , XcO, Xc2, and ir° . The results for the n° , although 
I ^~ ' off the experimental rate by 13%, is much closer than the usual expectations from a potential model. 
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I. CONSTRAINT TWO-BODY DIRAC EQUATIONS FOR QED AND QCD 

^ ' The Bethe-Salpeter equation (BSE) can be derived from QFT but its implementation invariably involves a non- 
. unique choice from a very large class of three dimensional truncations, almost all of which work well in the perturbative 
context. We describe one such approach, the Two-Body Dirac Equations (TBDE), based on Dirac's constraint 
dynamics that have been successfully applied to two-body bound state problems in QED [l[ and QCD Q and to two- 
body nucleon-nucleon scattering @ . We will describe various nonperturbative tests we have applied to distinguish it 
from the others @]. We then describe its most recent test in the treatment of two-photon decays of positronium and 
• quarkonium 0]. 

The TBDE can be derived from the BSE but they had their origins in classical relativistic mechanics [5|. For spin- 
zero particles one constructs two generalized mass shell constraints Hi = p 2 + mf + « 0; i = 1, 2 and guarantees 
their compatibility by requiring that the potentials satisfy a relativistic third-law condition $i = $2 = Pi, P2) = 

Oh, Q w and depend only on the interparticle separation component perpendicular to the total momentum, 
<D , 

x* 2± = (rr + P' l P")(xi-X2U P»=pZ+P% ; P"= P'VV^P 1 ; x 12 ± ■ P = 0. (1.1) 
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Thus, interactions depend on the invariant \/x±_ = \j x 2 + (x ■ P) 2 = r which reduces to the spatial separation 
in the cm frame, where the relative time cancels out covariantly. Calculating the difference of these compatible 
constraints, "Hi — H.2=— 2P ■ p ~ 0, gives a covariant elimination of the cm. relative energy, complimentary to the 
relative time restriction. The relative momentum p^— (s2P\ — £iP2)/ w IS given in terms of the total cm. energy 
w = E\ +£2 and the cm. constituent energies £1 — £2 = — ra^/w. It is canonically conjugate to x± in a covariant 
three-dimensional way, {x^_,p v } = T)j_. The other independent combination of the two constraints yields a dynamical 
equation (e2H\ + E\H-i)lw = p\ + — b 2 (u>) « whose quantization provides a covariant Schrodinger- like equation 
(p\ + & w )ip = b 2 (w)ip with the triangle function 

b 2 (w) = {w 4 - 2w 2 {m\ + m 2 2 ) + {m\ - m 2 2 ) 2 )/4w 2 = e 2 w - m 2 w (1.2) 

playing the role the eigenvalue. Its appearance signals exact two-body kinematics with effective particle motion 
displaying an Einstein relation between energy e w — {w 2 — m\ — m^jlw, mass m w = rnxm^jw, and momentum. 

For quantum mechanical spinning particles one has two Dirac equations (here given for minimal vector and scalar 
interactions) instead of two generalized Klein-Gordon-like mass shell constraints. 

Siip = 75i(7i ' (Pi - A*) + mi + Si)ip — 0; i = 1,2. (1.3) 

As in the spinless case they are compatible, [Si,S2]ip = 0, provided that supersymmetry is added to the conditions 
that apply in the two body spinless case. One finds that the constituent scalar Si = Si(S(r), A(r),p± , w, 71, 72) and 
vector = A^(A(r),p±, W, 71, 72) potentials are spin-dependent with the dynamics arising from single invariant 
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functions, one for each type of interaction. 

One of the advantages this covariant three dimensional equation has over the Bethc Salpeter equation is its simplicity. 
We obtain an invariant three dimensional Schrodinger-like equation 

{p 2 + a 2 ,p±, A(r), S(r))}V = & 2 (u# (1-4) 

in terms of a sixteen component wave function i/j = [?/>i, ip2, ^3: ipi]- Explicitly one finds that the equation has the 
following form exhibiting coupling between the upper-upper and lower-lower components in addition to the usual 
spin-dependent interactions symbolized by <&(<7i, 02,f>j_, A(r), S(r)), 

{p 2 + 2m w S + S 2 + 2e w A - A 2 + $ u }Vi + *i4^4 - b 2 {w)ijj = 0. (1.5) 

This equation and the one for the lower-lower component are quantum mechanically well defined, allowing nonpertu- 
bative numerical solutions. A second advantage is that in the case of lowest order QED where A(r) = ~a/r they have 
an analytic Sommerfeld-like solution for singlet positronium with spectral results that agree with standard perturba- 
tive results with the same spectral agreement holding numerically with triplet states This implies that they are 
less likely to produce spurious physics when applied to QCD as may occur in formalisms that do not have this type 
of nonperturbative QED agreement. Without successful nonperturbative QED tests how can a candidate two-body 
formalism be trusted for QCD spectral results? 

For QCD bound states we use a covariant version of the Adler-Piran static quark potential [|| Vap{t) = A(U(Ar) + 
Uq) (= A + S) obtained from an effective non-linear Maxwell equation embodying QCD. It analytically displays both 
asymptotic freedom Ai7(Ar << 1) ~ 1/rlnAr and linear quark confinement, AU(Ar > 2) ~ A 2 r. We obtain a 
covariant reinterpretation of their static model by replacing r with x± and apportioning the potential between our 
vector and scalar invariants A and S so that at short distance it is vector and long distance it is scalar Q. Once the 
underlying vector and scalar invariants and quark masses are fixed so are all the spectral predictions from Eq. (jl.4p . 

We obtain very good results for the entire meson spectrum from the light pion to the heavy upsilon states. The 
quality of our fit is about the same obtained by Godfry and Isgur Q, but with just 2 invariant functions instead 
of their 6 Q ■ The relativistic coupling structure of our equations are equally important for positronium and tt — p 
system. If we ignore the coupling to the lower-lower components ip^, we obtain poor hyperfine splitting results for 
the positronium system. Likewise we would obtain poor results for the hyperfine n— p splitting if we ignored that 
coupling (m,r ~ 850 MeV; m p ~ 1060 MeV) instead of the fully coupled results m v ~ 144 MeV; m p ~ 792 MeV @. 

The nonperturbative structures in our equations provide for chiral symmetry in that the pion (although not its 
excited states or p) behaves like a Goldstone boson. With the coupling structure, the pion mass tends to zero as 
the quark mass tends to zero, but without the coupling to lower-lower component this does not occur. However, 
the interaction structure of the TBDE when restricted to vector and scalar interactions does not give the functional 
dependence of on the quark mass of m 2 ~ m q behavior dictated from the nonconserved axial current generator 
of chiral symmetry. 



II. TWO-PHOTON DECAYS FROM TWO-BODY DIRAC EQUATION WAVE FUNCTIONS 

We take advantage of this close relation between the pion and positronium to treat two photon meson and positro- 
nium decays equally. The on-shell Feynman amplitude for electron-positron annihilation in the singlet state is 

M » s - (i^bf^^* - ^-Ziw-'ffl 7 ■ '"" + (1 * nut "' ip - > - (s+ w (2 - 1) 

The two-step change from above amplitude to that for positronium annihilation in which constituents are off shell 
is indicated below 

M aP ^ J d 3 P ^(p)M af3 = Mi So ^^ J d 3 ^(p)w(-p)r(p,k)u(p) 

= J d 3 P TP(p)TrT(p,k)u(p)v(-p) ^ J d 3 pTrT(p,k)V(p) (2.2) 

in which in the final step the 4x4 matrix wave function ^> replaces the outer product of the free u, v spinors. 
In Schr"odinger-like form the TBDE for singlet positronium becomes 



(-V 2 - 2e w a/r - a 2 /r 2 )ip(r) = b 2 {w)ij{r), 



(2.3) 
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with an exact solution wave function mildly singular at the origin. However, unlike standard approaches with 
rates proportional to ^(O)! 2 (which would give divergent results) the configuration space form of Eq. (|2.2p below 
smears the annihilation amplitude over a Compton wave length so that singularities are rendered harmless: 

- e^/^expHk ■ r)Tr { ^)b • ^m + ry • V) 7 ■ e 2+1 ■ e 2 {m + i, • V) 7 • ex P (-mr)/r}. (2.4) 

This gives rise to standard decay width formulae for positronium: 

T( 1 5 - 2 7 ) = ma 5 /2, 

T( 3 P -> 2 7 ) = 3ma 7 /256, T( 3 P 2 -> 2 7 ) = ma 7 /320, 

r( 3 p -» 2 7 )/r( 3 P 2 -» 27) = 15/4. (2.5) 
We recast the equations S,-^ = in terms of the mass and energy potentials 

Mi = rrii cosh L(S, A) + mi sinh £(5, A) 

E i = e l cosh£(A) -eiSinh^(A); G = exp£, (2.6) 

for which 

St = G&E, • p + EiPau + M a5 i - (i/2)GSj ■ 0(0$ + £/3 2 ) 7 5i 7 5 2 ; t, j = 1, 2. (2.7) 

For the pion and other 2- 7 meson decays we need to transform from 16 component wave function solutions of TBDE 
to their 4x4 matrix forms. The computation is simplified by first obtaining decoupled Schr"odinger like form of the 
TBDE through use of the plus/minus combinations <f>± = ipi ± ^4; x± = "02 ± "03 °f the original TBDE wave 
functions. A further spin-dependent scale transformation to 4 } ±tV± by the cm. transformation 

<t>± = exp(F + Ka- V rcr r r)ip ± ; x± = exp(F + K(T r v<T2-v)r)± (2.8) 

gives a still simpler decoupled form (with no momentum-dependent interactions), 

{p 2 + 2m w S(\r\) + S 2 {\r\) + 2s w A(\r\) - A 2 (\r\) + 9 D (S, A., VS, VA, V 2 S, V 2 A) 

+L • (<7i+<7 2 )$so(--) + L • (<ti-<t 2 )$sod(--) + iL-criX(T2$sco(~) 

+ cri-(T2®ss{--) + <ri-f<T 2 -f$ T (..) -I- cr r fcr2-fL • (<t 1 +<t 2 )$sot(--)}V'+ = b 2 (w)ip + . 

This four-component Schrodinger-like form has familiar spin-dependent interactions depending on the underlying 
invariants S and A that govern the scalar and vector interactions. Each of the four-component spinor wave functions 
is replaced by a 2x2 matrix wave functions 

<t>± = </>±ol + 4>± ■ X± = X±ol + X± ■ ° ( 2 -9) 

given in terms of singlet (scalar) {4>o± , Xo± > V'o± , Vo± ) an d triplet (vector) (4>±,x±,ip±,'n±) wave functions. 
The decay amplitude requires the full 4x4 matrix wave function 

^(r)=exp(F)[cosh.fir*(r) - sinh K S ■ f*(r)£ • f] 

\E r (r)=(V>+gi + ip-iq 2 + r/ + q +i]-q 3 ) /(2V2); q^ = 5. l;j q + ie ijk qk (2-10) 

Once we solve for four component spinor ip + we can use the two-body Dirac equation to solve for the remaining 12 

components, , if>_, T]o+,— W-i j with the final connection to the numerical meson wave functions from writing \l/(r) 

in terms of scalar Yj m and vector spherical harmonics Yj m ±,Xj TO , 

i'+ = 0+0 1 + ip+'tr 

0+o = {u%j/r)Y im ; V+ = (u+ +llj (r)/r)Y jm+ + ( U +_ nj (r)/r)Y jm _ + (uf +llj (r)/r)X jm -. (2.11) 
Of crucial importance is that the 16 component or 4x4 matrix form of TBD wave function must satisfy CM energy 
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dependent norm: 

' d 3 2# f (l + Aw 2 (3 1 f3 2 dA/dw 2 )ip} = I d 3 x^Cip = 1, (2.12) 



with the matrix A(r) =751752 — 71 • 72^]/2 displaying the core scalar and vector interactions. In 4x4 matrix form 
we have 

V'=exp(F)[coshif*(r) - sinhA'E ■ f*(r)£ • f] = /C*(r), (2.13) 
giving us the norm in terms ip and the Pauli reduction solution W(r). Thus in place of the Naive Norm (NN), 

J d 3 xTr¥ (r)^(r) = 1, 

we have the Two-Body Dirac Norm (TBDN), 

J d 3 xTr^&p = J d 3 xTr(K.y(r)) f OC*(r) = 1. (2.14) 

Using the (TBDN) norm has a significant effect on the decay amplitude and rate compared to that of the (NN) norm. 

For pseudoscalar meson decay, the triplet wave function is zero (ip + = 0) while for vector wave functions the singlet 
wave function (?A+o = 0) is zero. The final step of the decay formalism consists of substituting the resulting forms 
(I2.10p into Eq. (|2.4p for the matrix element and performing the numerical integration. 

Expt. TBDE(TBDN) TBDE(NN) 

Tr^ 1 ^- 0.135) 7.72±.04eV 8.73 eV 33.5 eV 

^(1%- 2.976) 7.4±1.0keV 6.20 keV 6.18 keV 

^(2% - 3.263) 1.3±0.6keV 3.36 keV 1.95 keV 

Xo(l 3 fb - 3.415) 2.6±0.65 keV 3.96 keV 3.34 keV 

X2(l 3 P2 - 3.556) 0.53±0.09 keV 0.743 keV 0.435 keV 

TABLE I: Meson 2 7 Decay Rates 

The results are for the ir° is very encouraging being off by only 13%. This is in sharp contrast to the 3 order 
of magnitude errors in most potential model applications. If we use the naive norm (NN) and amplitude then we 
obtain 33.5 eV. The results for the r\ c meson is about the same degree of accuracy. Work by Ackleh and Barnes Q 
and earlier ones by Haynes and Isgur [9} found it necessary to include results that appealed to effective field theories 
for the decay instead of the strictly microscopic approach we have taken here. 

In summary, the constraint TBDE formalism gives a covariant two-body wave equation that a) provides com- 
prehensive account for entire meson spectrum, b) is rigorously tested in QED, c) displays a remarkable connection 
between pion and singlet positronium that accounts for Goldstone boson-like behavior of pion, d) leads to a formalism 
that works well for light and heavy two-photon meson decays. 
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